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^ ; ABSTRACT 

. Context. Direct imaging of circumstellar disks at high angular resolution is mandatory to provide morphological information that 

^ ■ bring constraints on their properties, in particular the spatial distribution of dust. This challenging objective was for a long time in 

^ ' most cases, only within the realm of space telescopes from the visible to the infrared. New techniques combining observing strategy 

I and data processing now allow very high contrast imaging with 8-m class ground-based telescopes (10""^ to 10~^ at ~r') and comple- 

. ment space telescopes while improving angular resolution at near infrared wavelengths. 

■ Aims. We carried out a program at the VLT with NACO to image known debris disks with higher angular resolution in the near IR 
than ever before in order to study morphological properties and ultimately to detect signpost of planets. 

1—^ \ Methods. The observing method makes use of advanced techniques: Adaptive Optics, Coronagraphy and Differential Imaging, a com- 

ply , bination designed to directly image exoplanets with the upcoming generation of "planet finders" like GPI (Gemini Planet Imager) and 

I^T^ ■ SPHERE (Spectro-Polarimetric High contrast Exoplanet REsearch). Applied to extended objects like circumstellar disks, the method 

• ' is still successful but produces significant biases in terms of photometry and morphology. We developed a new model-matching pro- 

\ cedure to correct for these biases and hence to bring constraints on the morphology of debris disks. 

C^. Results. From our program, we present new images of the disk around the star HD 32297 obtained in the H (l.Gjim) and Ks {l.ljjLm) 

JL ■ bands with an unprecedented angular resolution (~65 mas). The images show an inclined thin disk detected at separations larger than 

^ ' 0.5 - 0.6". The modeling stage confirms a very high inclination (/ = 88°) and the presence of an inner cavity inside tq ^ 110 AU. We 

also found that the spine (line of maximum intensity along the midplane) of the disk is curved and we attributed this feature to a large 
anisotropic scattering factor (g ^ 0.5, valid for an non-edge on disk). 

Conclusions. Our modeling procedure is relevant to interpret images of circumstellar disks observed with Angular Diff'erential 
Imaging. It allows both to reduce the biases and to provide an estimation of disk parameters. 

' ^ ■ Key words. Stars: early-type - Stars: individual (HD 32297) - Techniques: image processing - Techniques: high angular resolution 

^ ■ 

m : 

n 1. Introduction these specific disks (warp, annulus, gap, off'set) can be ac- 

, , , , . ^„ , , counted for by the mass and orbital properties of these asso- 

. . Many stars have photometnc IR excess attributed to the pres- ^j^^^^ ^^^^^^ ^-^^^^ ^^^^^^^^ ^^j^^^ ^^^^^ ^-^^^ 

g ; ence of circumstellar dust. Among the circunistellar disk evo- where planets are not yet imaged. Many disks are thought to have 

^ lution, the debris disk phase is the period where most of the -^^^^^^^ ^^^-^-^^ attributed to dust sweeping by planets. From a 

^ pnmordial gas has been dissipated or accreted onto already ^j^^^^ structures are more or less 

■ formed giant planets. The dust content responsible for the IR ^^^^j identified. Therefore, a morphological analysis is manda- 
. ^ . excesses is thought to be the result of regular mutual collisions ^^^^ ^-^^ evolution, planetary formation and the likely 
k>H ■ between asteroid-like objects and/or evaporation or cometary presence of planets 

^ " bodies (Wyatt 2008; Krivov 2010). This activity presumably HD32297 is a young A5 star identified with an IR ex- 

_Cd. tiiggered by unseen planets is continuously replenishing the ^^^^ ^ ^ ^^^^^ ^^^^j^^^ ^ 

system with observable dust grains that, in turn, trace out the ^ j . ju- jii. oi j .1 /or^ricx • 

r 1 . . T-i . J . u J- . clmed circumstellar debris disk by Schneider et al. (2005) m 

presence or a planetary system. The recent detection by direct TVTTi^A/r/-.o ttot- .-.i- • .1 • -ui 

. . rt^. . .11 ji J • the near IR with NICMOS on HST. Observations m the visible 
imaging or planets m circumstellar disks around mam sequence 

A- type stars like j3 Pictoris (Lagrange et al. 2010) and HR 

8799 (Marois et al. 2008) or the undefined point source around 



(R band) were obtained by Kalas (2005) and reveal a circum- 
stellar nebulosity at distances larger than 560 AU indicating a 

Fomalhaut (Kalas et al. 2008; Janson et al. 2012) has clearly re- ^^^'^ scattered-light region extending to 1680 AU. The bright- 

, J . 1 • . w 1 . 1 \- ness asymmetries and a distorted morphology suggests the mter- 

vealed the suspected intimate relation between planet formation ^. vwi • . . n j- ttot- j ^ uf- ^ .^ ^ 

, , , , ^n.r^f\ or^r^o action With thc mtcrstellar mcdium. HST data obtamcd at 1 . liU m 

and debris disks structures (Ozernoy et al. 2000; Wyatt 2003; u.i- • ui Ar^r^ xtj j i o^iatt 

u . 1 ^nno\ A/r . • . / . • probc the inner region below -400 AU and as close as 34 AU 
Reche et al. 2008). Moreover, the most prominent structures in ^ , . ^ ^ i o j o^.x 
^ from the star (corresponding to respectively 3.3 and 0.3 ). 

Send offprint requests to: A. Boccaletti, e-mail: Schneider et al. (2005) derived a Position Angle (PA) of 47.6 ± 

anthony.boccalettiQobspm.fr 1.0° and an inclination of about 10° from edge-on. A surface 

* Based on observations collected at the European Southern brightness asymmetry is also observed in the South West (SW) 

Observatory, Chile, ESO (385.C-0476(B)) side of the disk which is about twice brighter than the North 
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Table 1. Log of observations summarizing the parameters: the filter used, the time of observation, the integration time of a single 
frame (DIT), the number of frames per data cube (NDIT), the number of cubes (Nexp), the seeing variation, the mean correlation 
time To, the amplitude of the parallactic angle variation and the observation time on target. 



Date 


Filter 


UT start/end 


DIT (s) 


NDIT 


Nexp 


seeinj 


?(") 


To (ms) 


Parall. amplitude (°) 


on target (s) 


2010.09.25 


H 


07:23:40/08:45:15 


2/3 


15/10 


68/13 


0.70- 


1.28 


1.7 


24.81 


1710 


2010.09.25 


Ks 


09:04:08 / 10:00:57 


2 


15 


57 


0.59- 


0.88 


2.0 


24.18 


2430 



East (NE) side at a projected angular separation of about 0.3- 
0.4". Observations at mid-IR wavelengths (11.7 and 18.3yum) 
by Moerchen et al. (2007) indicate a dust depletion closer than 
-70 AU suggestive of a inner cavity like the one found in 
HR4796 (Schneider et al. 1999) or HD 61005 (BuenzH et al. 
2010) but with a less steep decrease of the density. However, 
this radius is not well constrained in the mid IR as the equivalent 
resolution corresponds to 35 to 50 AU at the HD 32297 distance. 
Additional colors were obtained by Debes et al. (2009) at 1.60 
and 2.05yum again with NICMOS which confirms the disk ori- 
entation (PA = 46.3 ± 1.3° at NE and 49.0 ± 2.0° at SW). The 
surface brightness profiles confirm the asymmetries and reveal a 
slope break at 90 AU (NE) and 1 10 AU (SW). If similar to some 
other disks like j3 Pictoris (Augereau et al. 2001) and AU Mic 
(Augereau & Beust 2006) this break indicates the boundary of a 
planetesimals disk, it is marginally compatible with the value 
(70 AU) derived by Moerchen et al. (2007) but the measure- 
ments may not be directly comparable to that of Schneider et al. 
(2005). Also, it is shown that the disk midplane is curved at sep- 
arations larger than 1.5" (170 AU), possibly a result of the inter- 
action with the interstellar medium. This probably impacts the 
PA measurement. More recently, Mawet et al. (2009) achieved 
very high contrast imaging of the disk with a small (1.6 m) well 
corrected sub-aperture at the Palomar in Ks and confirmed the 
colors and brightness asymmetry, but the lower angular resolu- 
tion did not allow a detailed morphological study. 

At a distance of 1 13 ± 12 pc (Perryman et al. 1997), the disk 
of HD 32297 is of particular interest to probe the regions where 
planets would have formed in the 10- AU scale providing high 
angular resolution and high contrast are achieved simultane- 
ously. In this paper, we present new diff'raction limited images 
of the HD 32297 debris disk in the near IR (1.6 and 2.2 yum) 
using an 8-m class telescope with adaptive optics, therefore 
with a significant improvement with respect to former obser- 
vations (a factor of ~3 for the angular resolution at the same 
wavelengths). We note that similar observations obtained at the 
Keck telescope in the same near IR bands are presented by 
Esposito et al. (2012). Section 2 presents the results of the ob- 
servation. In section 3 we analyze the reliability of a particular 
structure observed in the image. Section 4 provides the Surface 
Brightness profiles of the disk before biases are taken into ac- 
count. The core of this paper is presented in the section 5 where 
we attempt to identify the parameters of the disk. In section 6 
the Surface Brightness corrected for the observational bias is de- 
rived. Finally, we present constraints on the presence of point 
sources within the disk. 



2. Observation and data reduction 

2.1. Observations 

The star HD 32297 (A5, V=8.13, H=7.62, K=7.59, near-IR 
values are from 2MASS) was observed with NACO, the 
near IR AO-assisted camera at the VLT (Rousset et al. 2003; 



Lenzenetal. 2003), on Sept. 25th, 2010 UT using a dedi- 
cated strategy to allow for high contrast imaging. We com- 
bined two modes off'ered in NACO: the four-quadrant phase 
mask coronagraphs (FQPMs, Rouan et al. 2000) to attenuate the 
on-axis star and the pupil tracking mode to allow for Angular 
Diff'erential Imaging (ADI, Marois et al. 2006). We used the 
broad-band filters H and Ks for which the FQPMs are optimized 
(Boccaletti et al. 2004) and the S27 camera providing a sampling 
of -0.027" per pixel. Coronagraphic images are obtained with a 
field stop of 8 X 8" and 13 x 13" in respectively H and Ks and 
an undersized pupil stop of 90% the full pupil diameter. Phase 
masks already provided interesting results on both extended 
objects (Riaud et al. 2006; Boccaletti et al. 2009; Mawet et al. 
2009) and point sources as companion to stars (Boccaletti et al. 
2008; Serabyn et al. 2010; Mawet et al. 2011). 

The log of observations is shown in Table 1. The integra- 
tion time of individual FQPM exposures (DIT) was chosen to fill 
about 80% of the dynamical range of the detector (including the 
attenuation by the coronagraph), and the duration of a template 
(DIT X NDIT) do not exceed 30s so as to correct for a system- 
atic drift correlated with the parallactic angle variations in pupil 
tracking mode. Therefore, the star is re-aligned with the corona- 
graphic mask every 30s when needed (the drift increases closer 
to the meridian). At the end of the sequence, a few sky frames 
were acquired around the target star. As for calibration purpose, 
we obtained out-of-mask images of the star itself with a neu- 
tral density filter (but without the undersized stop) to determine 
the point- spread function (PSF) which is intended to normalize 
coronagraphic images and hence to calculate surface brightness 
or contrast curves. Also, since the phase masks are engraved on a 
glass substrate, their transmission map is measured with a flat il- 
lumination provided by the internal lamp of NACO. Seeing con- 
ditions were good to medium: 0.72+0.06" in Ks and 0.97+0.12" 
in H. The AO correction provides a FWHM (Full Width at Half 
Maximum) of ~65mas in both filters while theoretical values for 
an 8-m class telescope at the longest wavelengths of the H and 
Ks bands are 46 and 59 mas respectively. Therefore, images are 
nearly diff'raction-limited in Ks as opposed to H. 

In pupil tracking, while the field rotates, the smearing of off'- 
axis objects is faster at the edge of the field. However, for a few 
seconds of integration per frame, the smearing is totally negligi- 
ble within our field of view. 

2.2. Data reduction and results 

The data were reduced with our own pipeline. First, we ap- 
plied bad pixel removal and Flat Field correction to every frame. 
Although quite faint here, we averaged sky frames and subtract 
a background map (including the dark) as it reduces the off'set. 
Given the large amount of coronagraphic images and to save 
computing time we binned the data cubes to reduce their size. 
We adopt a binning step of 5 frames (averaging) while keep- 
ing the ADI-induced smearing still negligible. A median sky 
frame is subtracted to every binned frames. Then, we rejected 
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Fig. 1. Final images of the disk around HD 32297 in the Ks (top) and H (bottom) bands for the cADI, rADI, and LOCI (left to right) 
algorithms (see section 2.2 for acronyms). The field of view of each sub-image is 3x3". The color bar indicates the contrast with 
respect to the star per pixel (not applicable to the display- saturated central residuals). 



the frames of worse quality with a flux criteria (open AO loop for 
instance) and obtained data cubes of 129 (75% of data) and 158 
(69% of data) frames for respectively the Ks and H band data, 
representing a total integration time of 1290 and 1706 seconds 
respectively. The frame recentering is a non linear problem, be- 
cause the star moves behind the coronagraph so the intensity and 
shape of the residual pattern vary significantly. It is therefore not 
possible to apply standard centroiding algorithm and the issues 
are diff'erent than in saturated imaging (Lagrange et al. 2012a). 
In particular, the FQPM creates four peaks of unequal intensi- 
ties which are difficult to exploit to provide an astrometric signal 
since phase aberrations are also quite large and perturb the coro- 
nagraphic image. We tested several procedures, focusing on the 
fitting of the halo. The best result is obtained when the individ- 
ual frames are thresholded to 10% of the image maximum and a 
Moff'at profile is fitted so that to give more weight to some pixels 
in the stellar halo with respect to those that are closer to the star. 
We achieved this way a sub-pixel accuracy for centering but we 
have not reliable mean to estimate a value. Finally, we used the 
background image of the FQPM alone to apply a correction of 
its own Flat Field to every re-centered frame. 

Following the nomenclature described in Lagrange et al. 
(2012a) we processed the data with several ADI algorithms 
(which diflfer in the way to build a reference frame): 

- classical ADI (cADI, Marois et al. 2006) which subtracts the 
median of the datacube to each frame. 



- radial ADI (rADI, Marois et al. 2006; Lafreniere et al. 2007) 
for which a set of reference frames to be subtracted is de- 
termined separately on each annulus. The algorithm is con- 
trolled by the separation criterion which defines in a deter- 
ministic way, together with the amount of field rotation, the 
time lapse between one frame and its reference frames. 

- Locally Optimized Combination of Images (LOCI, 
Lafreniere et al. 2007) is more sophisticated since it defines 
optimization zones (with a particular geometry) where the 
reference frame to be subtracted (in a zone smaller than the 
optimization one) is estimated as a linear combination of 
all frames where the field object has moved by the desired 
quantity (separation criterion Ns). Optimization zones are 
larger than the subtraction zones to limit the subtraction 
of a potential point source. Several parameters control the 
LOCI algorithm but the most relevant here are the separation 
criterion and the size of the optimization zone. 

After each algorithm, the frames are derotated to put the North 
vertical and a median is taken over to provide the final image. 

As circumstellar disks are extended objects, these techniques 
are not optimized since they favor the detection of unresolved 
sources and apply in diff'erent independent zones in the image 
breaking the continuous shape of a disk. Still, they can be very 
efficient to emphasize regions with strong gradients of similar 
spatial scales as point sources (rings for instance). In the past 
few years several impressive results were obtained using ADI 
and revealed the very detailed structures within some debris 
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disks (Buenzli et al. 2010; Thalmann et al. 201 1; Lagrange et al. 
2012a,b). Self-subtraction is an important issue and is caused by 
the contribution of the disk itself to the building of the reference 
image to be subtracted to the data cube. Because of the field rota- 
tion, the inner regions of a circumstellar disks are more aff'ected 
by self- subtraction in the ADI processing as they cover smaller 
angular sectors. Moreover, the self- subtraction of the disk in- 
creases as the algorithm ability to suppress speckles improves 
(LOCI being more aggressive than cADI for instance). This is 
particularly important in our data since the rotation amplitude is 
less than 25°. A detailed analysis of the self- subtraction arising 
when imaging disks with ADI will be presented in Milli et al. 
(submitted). 

ADI algorithms can produce a wealth of various images 
across which both the stellar residuals and the signal of the 
object will be modulated. The adjustment of the parameters 
relies on the balance between self- subtraction and the separa- 
tion from the star at which the disk is detected. This qualita- 
tive analysis leads us to set the separation criterion for both 
rADI and LOCI to 1.0 and l.SxFWHM for respectively the Ks 
and H bands, where FWHM=65 mas in both filters. The radial 
width is 2xFWHM for both the LOCI subtraction zones and 
rADI rings. As for LOCI, we adopt the geometry described in 
Lafreniere et al. (2007) where the zones are distributed centro- 
symmetrically. Subtraction zones have a radial to azimuth size 
factor of unity. The optimization zones start at the same radius 
than the subtraction zones but extend at larger radii with a to- 
tal area of Na = SOOxPSF cores. To build the reference image in 
each zone, all the available frames were used and the coefficients 
are calculated with a singular value decomposition. 

A first look at the resulting images shown in Fig. 1 indicates 
that 1/ no information is retrieved at separation closer than 0.5- 
0.6" since stellar residuals still dominate the image (for instance, 
the cADI images show strong waffle pattern as an imprint of the 
AO system), 2/ the H band data are of lower quality than in Ks 
(as expected from AO performance, see section 2.1, and seeing 
variations, see Table 1) but LOCI efficiently reduces the stellar 
residuals at similar radii in both bands, 3/ the disk has about the 
same thickness through all algorithms especially in the Ks band 
which favors a thin and highly inclined disk, 4/ the H and Ks 
morphologies are very similar. The structure of stellar residuals 
changes according to the algorithm, and diff'ers notably from the 
well-known speckled pattern, especially owing to the azimuthal 
averaging caused by the de-rotation of the frames. A closer in- 
spection shows that although the disk is quite linear it features a 
small deviation from the midplane which is more emphasized 
in LOCI images to the NE side (especially in Ks but also in 
H). Figure 2 shows the position of the midplane (red line) with 
respect to the spine of the disk. The deviation is evidenced in 
the middle of this image and spans several resolution elements 
(from -0.6" to -1.0"). The main objective of the paper is to find 
a disk morphological model which can account for this devia- 
tion. In addition, we present signal to noise ratio maps in Fig. 
3. The noise level was estimated azimuthally on LOCI images 
(both H and Ks) which are smoothed with a gaussian (FWHM 
= 65ms) so as to derive the SNR per resolution element. We 
achieve detection levels of about 10 to 20cr along the disk mid- 
plane. We warn that the noise estimation can be biased especially 
in the speckle-dominated regime close to the star as explained in 
Marois et al. (2008). In the case of Fig. 3 this occurs at distances 
closer than 0.5 ~ 0.6" although a few speckles still appear fur- 
ther out as significant pattern in the H band SNR map. Finally, 
we note that the deviation from the midplane is also visible in 
the SNR maps. 



Fig. 2. Upper left 1x1" quadrant of the Ks LOCI image which 
emphasizes the deviation with respect to the midplane (red line, 
PA = 47.4°). The star is at the lower-right comer. The color 
bar indicates the contrast with respect to the star per pixel (not 
applicable to the display- saturated central residuals). 
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Fig. 3. SNR map obtained from LOCI images in Ks (top) and 
H (bottom) for SNR > 5cr per resolution element. The field of 
view is 4"xl.5" and the disk is derotated (PA = 47.4°) to put the 
disk horizontal. The outer (> 1.5") and inner (< 0.6") regions 
are excluded. The color bar indicates the SNR per resolution 
elements. 



3. Morphological analysis 

3.1. Position Angle determination 

First, we measured the PA of the disk to check its consistency 
with the values reported previously by Schneider et al. (2005) 
and Debes et al. (2009), and with the goal to improve the accu- 
racy owing to the higher angular resolution and the benefit of the 
ADI processing. For that, we followed the approach presented 
in Boccaletti et al. (2009) and more recently in Lagrange et al. 
(2012a) for the case of JS Pictoris, as it is nearly edge-on as the 
one around HD 32297. The image is de-rotated by an angle cor- 
responding to a rough estimation of the PA (guess) so that the 
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disk is about horizontal with respect to pixels. Vertical fitting 
with a gaussian profile along the columns yields for each sepa- 
ration the position of the spine (trace of the maximum surface 
brightness along the midplane) at a sub-pixel precision. Then, 
from the position of the spine as a function of radius (smoothed 
over 4 pixels ^2 FWHM), we estimated the slope of the spine 
in four regions carefully chosen to reduce the influence of stellar 
residuals and background noise. The measurements were made 
in the regions 25-40, 25-50, 30-40 and 30-50 pixels (correspond- 
ing to a range of 0.675 - 1.35") on each side of the disk together 
(assuming the disk crosses the star). 

The process is repeated for several angles in a range of 
2° sampled at 1/100° around the guessed PA. Finally, the spine 
slope as a function of this angle shows a steep minimum which 
defines unambiguously (at a precision of 0.2°) the PA of the disk. 
We made the same analysis for all the images in Fig. 1 and aver- 
aged the values. The error bar is calculated from the dispersion 
of the values to which we added quadratically the Celestial North 
uncertainty (0.2° according to Lagrange et al. 2012a). We found 
PA = 47 A ± 0.3° in the Ks band and PA = 47.0 ± 0.3° in the 
H band which agree within error bars but we adopt the former 
value since Ks data are of much better quality. These numbers 
compare adequately to previous reports but achieve better pre- 
cision. However, we did not observe a PA variation of about 3° 
between the two sides of the disk unlike Debes et al. (2009) but 
we perform our measurement at closer separations than in HST 
so the studies are not directly comparable. 

3.2. Reliability of the spine deviation from midplane 

Following the observation reported in the previous section, our 
main objective was to test if the deviation of the midplane seen 
in Fig. 1 and Fig. 2 could be a bias produced by the ADI al- 
gorithms. Therefore, to overcome this issue and to transform 
observed measurements into physical ones requires calibration 
through models. In view of our images, unlike HST ones, we 
are not sensitive to the regions where the disk is posited to in- 
teract with the ISM (because of a lower sensitivity). Hence, we 
assumed a simple geometrical model of a ring-like dust disk. 
We used the GRaTer code (Augereau et al. 1999; Lebreton et al. 
2012) to calculate synthetic scattered light images of optically 
thin disks, assuming axial symmetry about their rotation axis, 
and inject such models in the real data. In this study, we assume 
a dust number density that writes (in cylindrical coordinates) 



-2ain I ^ \ -2aoi 



(( 



exp 



(1) 



where and //q are the midplane number density and the disk 
scale height at the reference distance to the star r = ro, respec- 
tively. The slopes Qf/„ and aout are usually chosen such that the 
density is rising with increasing distance up to ro and is decreas- 
ing beyond {a in > and aout ^ 0)- This is aimed to mimic a 
dust belt with a peak density position close to ro, a geometry ap- 
propriate for cold debris disks in isolation (no interaction with 
ISM). A ray-traced image is then calculated assuming a disk 
inclination / from pole-on viewing, and anisotropic scattering 
modeled by an Henyey & Greenstein (1941) phase function with 
asymmetry parameter g between - 1 (pure backward scattering) 
and 1 (forward scattering). To limit the number of free parame- 
ters, we assumed a gaussian vertical profile (y = 2) with a scale 
height proportional to the distance to the star (J3 = 1). At this 




radius in arcsec 



Fig. 4. Width of the observed disk (black line) in Ks for cADI 
(top), rADI (middle) and LOCI (bottom) images, compared to 
the width of the model for / = 80° (blue line), / = 85° (magenta 
line) and / = 88° (red line). The NE side is at the left. 



stage, we generated a limited number of models with some as- 
sumptions on the remaining free parameters taken from the liter- 
ature (Schneider et al. 2005; Moerchen et al. 2007; Debes et al. 
2009): g = (isotropic scattering), ro = 70 AU, aout = -4, 
Ho = 2AU, / = [80,85,88]°, atn = [0,2]. The vertical height 
Hq is consistent with the natural thickening of debris disks as 
defined by Thebault (2009), although chosen close to the lower 
boundary of the proposed range in that paper. Hence in our case, 
the aspect ratio is assumed constant: H/r ^ 0.029. We note that 
the simplified version of GRaTer used in this study does not han- 
dle perfect edge-on viewing i=90°. 

Images of the models are first convolved with the observed 
PSF (obtained with the ND filter but no stop). While the field 
rotates, the orientation of the disk changes with respect to some 
PSF features (spiders and static phase aberrations). So rigor- 
ously, a model must be first rotated at a given orientation and 
then convolved with the PSF, and the convolution is repeated 
for each parallactic angle in the data cube. We checked however 
that our conclusions are not aff'ected if the convolution is ap- 
plied once and the convolved model is rotated. Likely, the AO- 
correction quality is not good enough in this particular case to 
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Fig. 5. Spine of the disk in Ks for the observed disk (bluish 
Hnes), and the / = 88° model (reddish lines) as measured in 
cADI, rADI and LOCI. The elevation from the midplane is mea- 
sured for PA = 47.4°. Arrows indicate where the deviation to 
the midplane is observed. The NE side is at the left. 



make a difference. We adopt this faster solution to save comput- 
ing time. Model disks were inserted at 90° from the observed 
disk so there is no possible overlap in between. A similar strat- 
egy was adopted by Thalmann et al. (2011). We used the same 
approach as described above for measuring the PA. The fitting 
of the observed and model disk vertical profiles versus radius 
yields the disk spine position, vertical thickness and intensity. 
Vertical thickness appears to be a good criterion to disentangle 
between several disk inclinations (/ parameter). Figure 4 plots 
the vertical thickness (FWHM of a gaussian) measured on the 
observed disk (black line) and the model disks (colored lines) 
for the three inclinations and the three algorithms. Qualitatively, 
we found that the best matching with the models is achieved for 
an inclination of / = 88°. This will be confirmed in the next sec- 
tion with a more precise model fitting but already seems diff'erent 
than the values of Schneider et al. (2005) where 77° < / < 82°. 
Also, we note that the shape of the FWHM profiles as a function 
of radius is quite well reproduced with the model, providing sup- 
port to a constant disk aspect ratio. This suggests that the disk is 
opened beyond about 1". Finally, we measure for this particu- 
lar inclination of 88°, a minimal observed vertical thickness for 
the observed and the model disks of about 0.07" at a distance 
of about 0.85" (96 AU) in cADI images. This is consistent with 
the value of Hq used in the model given the eff'ect due to the 
inclination and the convolution with the PSF. 

Figure 5 shows for the most likely inclination (/ = 88°) the 
disk spine position for cADI, rADI and LOCI images in Ks and 
for both the observed and model disks. The deviation to the mid- 
plane is evidenced especially in the NE side for Ks images be- 
tween 0.6 and 1" and reaches at maximum 0.025" (~1 pixel). 
On average the deviation has an amplitude of about 0.015" and 
spans across 15 pixels (which is about 6 resolution elements). 
The noise level per resolution element is equivalent to about 
0.005". Therefore, the deviation from the midplane is significant 
at the ~7.5-sigma level. A low-significance symmetrical devia- 
tion is barely seen in the SW side too with the same amplitude. 
In H-band images, the deviation is more difficult to detect. It is 
still visible in the NE side in rADI and LOCI images only. There 
is also a possible deviation at larger radii (> 1.2") on both sides 
and towards the North- West which may comply with the change 
of PA noted by Debes et al. (2009) but our data lack of sensi- 
tivity compared to HST in these regions for a confirmation. We 
checked that other models mostly lead to the same conclusion. 



Imaging of HD 32297 

Table 2. Slopes (as fitted to power laws) of the different com- 
ponents in the surface brightness of cADI images. 





Hband 


Kband 




NE 


SW 


NE 


SW 


inner slopes 


1.9 


1.4 


-0.6 


0.0 


outer slopes 


-3.4 


-3.2 


-4.6 


-3.9 



In addition, we found that the deviation in the NE side at 0.6-1" 
is resilient to a variation of the PA within the error bar (0.3°). 

Finally, as the LOCI performs locally the injection of model 
disks at 90° from the observed disk, the procedure might be in- 
appropriate to calibrate the bias of the algorithm. Instead, in the 
LOCI images of Fig. 1, we calculated the coefficients of the lin- 
ear combinations for each zone and apply them to a model disk 
at the same PA than the observed disk (but in an empty data 
cube) and still observe a nearly ffat midplane of the model disk 
in the 0.6-1" region. To summarize, an inclined model disk with 
a symmetrical vertical intensity distribution with respect to the 
midplane (g = 0, i < 90°) produces the same symmetry in the 
ADI processed images. Therefore, the shape of the midplane ob- 
served in Figs. 1, 2 and 5 is definitely not due to a bias of the 
ADI process but is real. Interestingly and as an additional evi- 
dence, Esposito et al. (2012) images from Keck reveal (at least 
qualitatively) the same deviation from the midplane at the same 
separations. 

4. Photometric analysis 

The photometric analysis of a disk observed in ADI can be prob- 
lematic especially when the field rotation is small as in our case. 
The ADI effect can be understood in the resulting image as a 
transmission map projected on the sky in which the pattern de- 
pends on the disk brightness distribution. It is for instance very 
interesting to note in Lagrange et al. (2012a) that LOCI produces 
on the disk of JS Pic a very strong attenuation, just because it is 
very well resolved in the vertical dimension. The actual bright- 
ness distribution of the disk is therefore needed to properly ac- 
count for this attenuation and derive reliable measures. This is 
the motivation for the modeling work in the next section. The 
Surface Brightness profiles presented in this section cannot be 
considered as the true photometry but are providing us with a 
starting point. 

As a first step, we performed the Surface Brightness (SB) 
measurement on cADI images as cADI produces less self- 
subtraction than other algorithms and hence less photomet- 
ric/astrometric biases. The cADI image of the disk is smoothed 
over 4 pixels, and isolated in a narrow slit of 0.25" in width 
aligned on the disk PA. From that, we extracted the mean in- 
tensity profile of the disk azimuthally averaged in angular sec- 
tors as well as the intensity dispersion and plot the results in 
Fig. 6. Doing so and given the edge-on geometry, the intensity 
mean and dispersion are hence calculated in pixels located at 
same physical distances to the star (same ffux received) but dif- 
ferent elevations from the midplane. In addition, we did the same 
measurement out of the slit to derive the noise level. There is a 
clear break in the SB profiles near 1" (about 1 10 AU) between a 
nearly ffat (or even positive) proffle inwards and a negative slope 
outwards. The slopes given in Table 2 are measured at speciffc 
positions that we consider to be the most relevant: 0.6-0.9" for 
the inner slopes and 1.2-1.8" for the outer slopes. Clearly, the 
inner slopes have lower signiffcance because 1/ they are mea- 
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Angular distance in arcsec 



17 
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Fig. 6. Surface Brightness expressed in mag/arcsec^ as measured 
in the H band (top) and in the Ks band (bottom). Red and blue 
Hnes stand for the NE and SW sides respectively. Dotted lines 
show the disk intensity dispersion at 1 cr within the slit window 
of 0.25". The black line indicates the noise level measured out 
of the slit window. Dashed lines correspond to the fitted power 
law profile next to which the slopes are indicated. 



sured across only 0.3", equivalent to 11 pixels (4.5 resolution 
elements), and 2/ this region has higher speckle noise (as indi- 
cated by the noise level in Fig. 6). Indeed, the error bar on the SB 
slopes can be large considering the intensity dispersion plotted 
in Fig. 6. To be conservative we adopt an uncertainty of ±0.5 for 
the index of the slopes. The inner slope is steeper in H than in Ks 
but this is not reliable given the noise level at such separations 
as seen in Fig. 1. The outer slopes are very similar for the two 
sides of the disk (within error bars), but noticeably steeper in 
Ks than in H as already noted by Debes et al. (2009). The disk 
also appears dimmer in H than in Ks which might be consis- 
tent again with HST data. As in the case of p Pictoris we may 
expect that the disk of HD 32297 contains a planetesimals belt 
at radii smaller than ~ 110 AU which produces dust particles 
blown away from the system by the radiation pressure of the star 
(Augereau et al. 2001). Outer slopes close to -4 are compatible 
with such an assumption. 



5. Disk modeling 

5.1. Parameter space 

We now explore if the morphology of the disk can be ac- 
counted for by a model of scattered light images. Using GRaTer 
(Augereau et al. 1999), we generated a grid of 576 models of 
which the parameters span relevant values (based on previous 
guess) for the disk of HD 32297 as given below: 

- inclination (/): 82, 85, 88, 89.5° 

- peak density position (ro): 80, 90, 100, 110, 120 AU 

- power law index of the outer density {a out)'- -4, -5, -6 

- power law index of the inner density (a/«): 2, 5, 10 

- anisotropic scattering factor {g)\ 0, 0.25, 0.50, 0.75 

We considered only positive values of g since we suspect 
that the deviation from the midplane is in fact due to anisotropic 
scattering towards the NW, but we cannot make any distinction 
with negative values and an inclination larger than 90°. In addi- 
tion to all possible combinations of disk parameters we have two 
filters and three diff'erent images (cADI, rADI, LOCI) for each, 
which totals 3456 prescriptions. Based on the analysis in sec- 
tion 3, we rejected an inclination of 80° and ro = 70 AU as they 
are clearly not compatible with our data. Similarly to section 3, 
the disk models were convolved with the PSF, injected at 90° 
from the disk PA and normalized with the observed disk in the 
outer part of the disk (40-50 pixels for Ks and 30-40 pixels for 
H). However, the intensity normalization is not extremely accu- 
rate for two reasons: given the possible range of parameters the 
fake disk morphology can be diff'erent than that of the observed 
disk, and also the ADI algorithms, especially LOCI, produce a 
self- subtraction that depends on the morphology. Therefore, in 
our eff'ort to match the observed disk to models we proceeded 
in several steps. We first started with a global minimization of 
images pixel to pixel and then we focused on minimizing partic- 
ular measurements like the slopes in the Surface Brightness, and 
the spine of the disk. At the end, we are looking for a set of disk 
parameters [/, ro, aout^ (^in, g] that matches as far as possible the 
observed disk. 



5.2. Global minimization 

As for the global minimization we isolated the relevant part of 
the disk with a slit centered on the disk PA respectively PA-h90° 
for the model, 0.5" in width for Ks respectively 0.3" in H and 
restrained the radius to the range 0.65-2" for H and Ks. These 
numbers were carefully chosen to retain only the pixels contain- 
ing flux from the observed disk and hence optimize the mini- 
mization. Then for each model we used a standard metric as 
follow: 



2 _ (.hiskiU j) - ImodeliU j)) 



ij 



cr- 



disk 



(iJ) 



(2) 



where Idisk and Imodei are the images of the observed disk 
and that of the model for the retained pixels (/, j), and cr^.^^ is 
derived from the azimuthal dispersion of Idisk to account for the 
errors. As we run the minimization in both sides NE/SW of the 
disk simultaneously, it implicitly assumes that the disk is sym- 
metrical about the minor axis. The;^^ minimum usually gives the 
adequate model that best matches the data, but to prevent degen- 
eracies between parameters of the model we selected instead the 
10 lowest minima corresponding to 10 sets of model parameters 
per ADI algorithm. Combining the three algorithms amounts to 
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Table 3. Best model parameters with the global minimization. Table 4. Best model parameters with the Surface Brightness 
Significance of values are indicated in parenthesis. minimization. Significance of values are indicated in parenthe- 

sis. 



filter / (°) Oout rp (AU) g 

H 88 (0.70) -6(0.43) 10(0.47) 110(0.47) 0.50(0.37) 
Ks 88(0.97) -5 (0.5) 2(0.37) 110,120(0.5) 0.50(0.90) 



30 sets of parameters. It would not be appropriate to derive av- 
eraged values nor error bars since our grid of models is discrete 
and cannot be interpolated. Instead, we searched independently 
for each parameter among these sets which value has the highest 
occurrence as a measure of the significance (value with the high- 
est occurrence/30). This quantity is a number between and 1 
and does not convey an assumption of the error bar. The results 
are given in Table 3 for both H and Ks. In the particular case of ro 
we found that both 110 and 120 AU were equally probable in the 
Ks image. The significance (indicated in parenthesis in Table 3) 
has to be larger than the inverse of the number of allowed values 
for each parameter. Some parameters are very well constrained 
like / and g (at least in Ks) supporting a quite large inclination 
of 88°. The value for ro is in very good agreement with the slope 
break in the Surface Brightness measured by Debes et al. (2009) 
as in our measurement (1 10 AU on average). Density profiles are 
more difficult to constrain especially Qf/„ where both H and Ks 
are not consistent but we may expect chromatic dependence for 
this parameters or NE/SW diff'erences. Finally, the global mini- 
mization clearly points to a noticeable anisotropic scattering of 
g = 0.5. Minimization in H and Ks leads to the same values for 
/, ro and g although the significance is always lower in H as a 
result of lower data quality. 

5.3. Minimization of Surface Brightness slopes and spine 

The global minimization approach is useful to get a global pic- 
ture of the disk because it takes into account all the available 
information in all the relevant pixels and try to simultaneously 
adjust all the parameters of the model, assuming the disk is sym- 
metrical with respect to the minor axis. But since the intensity 
scaling of the model with respect to the observed disk can be 
problematic (as explained in the previous sub- section) we tested 
alternative procedures. To refine the analysis, as a second step, 
we carry out the minimization with other criteria that are less 
dependent on pixel to pixel variations. Here, we set the inclina- 
tion to / = 88° as a clear outcome of the global minimization 
analysis. 

Starting with the Surface Brightness we proceed similarly to 
section 4 to extract this quantity on both the observed disk and 
the model disks. We do not minimize the full Surface Brightness 
profile but just the slopes in the 4 areas of interest inward and 
outward of 1" and independently in the NE and SW sides. We 
used the same regions defined in section 4 to minimize aout and 
ain, so respectively 1.2-1.8" and 0.6-0.9". Again, the minimiza- 
tion provides 30 sets of parameters for each area and we searched 
for the highest occurrence. The results are shown in Table 4. The 
analysis of aout shows a steeper index in the NE but differs for 
Ks and H but remains in agreement with the SB profiles of the 
observed disk plotted in Fig. 6. The inner density indices are dis- 
crepant with the global minimization, although with larger sig- 
nificance and rather converge to an intermediate value of Qf/„ = 5. 
The parameter ro is more subject to dispersion. We concluded 
that these values are not representative since the range of ro falls 



filter 


area 






ro (AU) 


g 


H 


NE in 
SW in 
NE out 
SW out 


-5 (0.50) 
-4 (0.57) 


5 (0.50) 
5 (0.50) 


120 (0.57) 
120 (0.47) 
120 (0.57) 
100 (0.40) 


0.25 (0.63) 
0.25 (0.50) 
0.75 (0.33) 
0.25 (0.30) 


Ks 


NE in 
SW in 
NE out 
SW out 


-6 (0.70) 
-5 (0.47) 


5 (0.47) 
2 (0.40) 


100 (0.27) 
90 (0.47) 
120 (0.30) 
100 (0.37) 


0.75 (0.30) 
0.50 (0.47) 
0.75 (0.73) 
0.00 (0.47) 



unequally between the two regions where the minimization is 
performed. However, in this particular case, it is informative to 
use a mean instead of the highest occurrence since the allowed 
values for ro are better sampled than for the other parameter. This 
mean is much closer to 110 AU in concordance with the global 
minimization. Furthermore, the minimization of SB slopes is not 
appropriate to constrain the anisotropic scattering factor since 
the dispersion of solutions is too large (0.00<g<0.75). 

One last measurement was carried out to further constrain 
the anisotropic scattering factor. We applied the procedure de- 
scribed in section 3 to evaluate the spine of the disk as we con- 
sider it can provide a better estimation for the anisotropic scat- 
tering factor although potentially correlated to other parameters 
like the densities. For the same reasons as above, we expect these 
geometrical measurements to be less dependent on the intensity 
normalization of the fake disks. Results are reported in Table 5 
where we considered the NE and SW sides separately. Here, we 
did not expect significant results for agut and a/„ which are re- 
lated to the disk intensity while we measure just coordinates for 
the spine. This is confirmed in Table 5 where no particular value 
emerges. The radius ro shows a difference in between the NE 
(1 10-120 AU) and SW (80-90 AU) sides (although with a low 
significance in Ks) which actually is in agreement with the im- 
ages where the NE side appear slightly more elongated than the 
SW side. Debes et al. (2009) also noted such a NE/SW asymme- 
try but in the opposite direction. However, they used directly the 
indication from the SB without modeling. Orbital motion of a 
hypothetical structure at such a separation is totally inconsistent 
between the two epochs. As the disk is seen almost edge-on, this 
asymmetry could be related to a stellocentric offset of the cen- 
tral cavity as in many other debris disks (Boccaletti et al. 2003; 
Kalas et al. 2005; Thalmann et al. 2011; Lagrange et al. 2012b) 
and then interpreted as a signpost of planet which is sculpt- 
ing the disk. The relation between this potential offset and the 
anisotropic scattering remains to be explored. As for the factor 
g, the Ks band favors one single value of 0.25, while it is larger 
in H. We note that the largest value of the anisotropic scattering 
factor (g = 0.75) is actually unrealistic according to our images 
since it would make the disk brighter than the stellar residuals in 
the zone r < 0.65". As we exclude this region dominated by the 
starlight in the minimization process, the relevant pixels where 
large values of g contributes are missed and hence some large 
values may come out. It is surprising that the spine minimization 
yields a different value for g than the global minimization (0.25 
instead of 0.50) while it remains symmetrical. Nevertheless, the 
global minimization converges to g = 0.50 with a much higher 
significance at least in Ks. 
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Fig. 7. Images of the best disk model (PA = 137.4°) in Ks compared to the observed disk (PA = 47.4°) with cADI, rADI, and LOCI 
algorithms (left to right). The field of view of each sub-image is 3x3". Color bars indicate the contrast with respect to the star (not 
applicable to the display- saturated central residuals). 



ro and g are derived with quite good significance, at least for the 
global and spine minimization. For some parameters, it is rea- 
sonable to assume a higher weight from the images in Ks which 
are of much better quality than in the H band. Considering an 
axisymmetrical disk, we retain the values of the global mini- 
mization ro = 1 10 AU and g = 0.50 as the optimal numbers 
to describe the disk morphology. More difficult to constrain are 
the indices of the density profiles. First of all, inward from ro 
the number of pixels used to fit a slope is very much reduced. 
The 0.6-0.9" region used corresponds to about 11 pixels only 
in radius. Second, the noise introduced by the residual, corre- 
sponding to un-attenuated starlight produces bumpiness in the 
Surface Brightness which can easily aff'ect the measurement of 
the slope exponent by 1, while it is the level of precision we 
are looking for as imposed by the model parameter sampling. 
Combining the global and Surface Brightness minimizations we 
choose aout = -5 in Ks. In H, the values of aout are discrepant. 
The SB minimization, in principle more accurate to measure the 
SB slopes, points to a lower value. In order to be consistent with 
the slopes shown in Fig. 6, we adopt a out = -5 in H. This value 
is consistent with the predictions of Lecavelier des Etangs et al. 
(1996) for a circular ring of planetesimals producing dust with 
collisions and expelled by the radiation pressure. It is also sim- 
ilar to the slopes measured in the debris disk around J3 Pic 
(Golimowski et al. 2006). Finally, the values of a/„ are in prac- 
tice clearly more difficult to measure but appear systematically 
larger in H than in Ks. To reproduce such a behavior we select 
a in = 5 in H and a in = 2 in Ks. 

Table 6 gives the final choice of parameters combining the 
results presented above in Tables 3, 4 and 5. As for illustration, 
a model disk with the parameters for the Ks band is shown in 
Fig. 7 at 90° from the observed disk. In both filters, the parame- 
ters that describe the morphology of the disk as projected on the 
sky (/, ro and g) are strictly identical. This is definitely a good 
indication that the minimization process is performing well. As 
an additional demonstration. Fig. 8 shows for the maximum in- 
clination allowed by HST measurement (Schneider et al. 2005) 
that the model would be discrepant with our disk image. If the 
inclination were / = 82° with no anisotropy, we should see either 
the inner cavity or the a thicker disk depending on a in. With an 
anisotropic scattering factor identical to our estimate (g = 0.5), 
the deviation with respect to the midplane would be larger that 
what we measured. The parameters that rather depend on inten- 



Table 5. Best model parameters with the spine minimization. 
Significance of values are indicated in parenthesis. 



filter 


area 


0!out 




ro (AU) g 


H 


NE 
SW 


- 4 (0.37) 

- 6 (0.73) 


10 (0.60) 
5 (0.37) 


120 (0.53) 0.5 (0.67) 
80 (0.43) 0.75 (0.70) 


Ks 


NE 
SW 


- 6 (0.40) 

- 4 (0.53) 


5 (0.47) 
2 (0.57) 


110(0.23) 0.25 (0.53) 
90 (0.33) 0.25 (0.67) 




Fig. 8. Images of the disk model (PA = 137.4°) in Ks for / = 82° 
and g = 0, 0.5 (respectively left and right) compared to the ob- 
served disk (PA = 47.4°) with rADI. The field of view of each 
sub-image is 3x3". Color bars indicate the contrast with respect 
to the star (not applicable to the display- saturated central resid- 
uals). 



5.4. Discussion 

The determination of the model parameters is not perfect in the 
case of our data given the signal to noise ratio of the images 
while we also expect some degeneracies in between some pa- 
rameters. The disk of HD 32297 being close to edge-on, the 
inclination is definitely the easiest parameter to measure. We 
found an optimal value / = 88°, larger than the original mea- 
sure with HST which makes the disk more inclined than pre- 
viously thought. This difference of inclination is possibly ex- 
plained and hence actually supports the scenario of interaction 
with the ISM (Debes et al. 2009) which has noticeable effect in 
HST data because of a better sensitivity to faint sources Then, 
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Table 6. Best model parameters retained. 



filter / (°) 


^out 




ro {AU) 


g 


H 88 


-5 


5 


110 


0.50 


Ks 88 


-5 


2 


110 


0.50 



sities {aout^ o:in) are on the contrary subject to questioning. The 
difference observed between the two filters can be due to either 
uncertainties or chromatic dependence. The accuracy of the data 
does not allow to firmly decide and these power low indices gen- 
erate SB slopes which agree within error bars with those mea- 
sured in Fig. 6. However, we note that the depletion of the central 
cavity is not as steep as other cases like HD 61005 (Buenzli et al. 
2010) or HR 4796 (Schneider et al. 1999; Lagrange et al. 2012b) 
as already mentioned by Moerchen et al. (2007) from mid IR im- 
ages. 



6. Surface Brightness of the best models 

The Surface Brightness profiles presented in Fig. 6 are impacted 
by the ADI process which produces a radial dependent self- 
subtraction of the disk. In particular, the parts of the disk that 
are closer to the star undergo a larger subtraction. 

However, since the best model has been identified we can 
attempt to recover the most likely SB profile corrected from the 
self- subtraction by measuring it directly on the model once nor- 
malized adequately. For that, we first compare the cADI image 
of the observed disk and the cADI image of the best model to 
refine the intensity ratio in between. We used this ratio to in- 
ject a fake disk in an empty data cube to mimic the field rota- 
tion. Instead of processing the data with one of the ADI algo- 
rithms we simply derotated the frames and sum them up. Then, 
we applied the same measurement as in section 4. The resulting 
SB profiles corresponding to the aforementioned parameters are 
shown in Fig. 9. The SB profiles derived from the models are 
significantly diflTerent than those measured in cADI images. ADI 
has an obvious ejfect on the slopes of the SB as also noted in 
Lagrange et al (2012b). The net result of the ADI with respect 
to the model, in our case, is to change the intensity variations 
along the midplane sometimes reinforcing the intensity gradi- 
ent. Starting from a negative inner slope in the models (Fig. 9) it 
becomes nearly flat (in Ks) or even positive (in H) in the ADI im- 
ages (Fig. 6). Similarly, the outer slope exponents are increased 
by ~ 1 . We were not able to reproduce the steepness of the inner 
profile in H but the negative areas in the images (Fig. 1) suggest 
a stronger eflfect of the reference frame subtraction than in Ks 
^ . Finally, the profiles for the models are brighter as expected 
since the ADI self- subtraction is accounted for. The cADI im- 
age in Ks is about 0.3 mag/arcsec^ fainter than the model at 1" 
(corresponding to ro). In H, the attenuation can be as large as ~1 
mag/arcsec^ at the same radius. 



^ In the H band, the decline close to the axis of the SB as observed 
in Fig. 6 is caused by the negative pattern that appear in Fig. 1 (bottom 
left) at a diff'erent PA than the disk but still aff'ects the disk. It is due 
to a poor estimation of the reference image essentially because of the 
diffraction pattern decorrelation over time. We suspect that this butterfly 
pattern (positive and negative) is due to the secondary spiders. The angle 
is very close to the 109° that separates two VLT spiders. This residual 
pattern is seen in several other NACO images. 



13 F ' ' ' ' ^ 

I NE side 

-1 4 L SW side 




Angular distance in arcsec 

13 F ' ' ' ' \ 

I NE side 




Angular distance in arcsec 



Fig. 9. Surface Brightness expressed in mag/arcsec^ as measured 
for the model with parameters / = 88°, ro = 110A[/, a out = -5, 
a in = 5, g = 0.50 (top), corresponding to the H band and for 
the parameters / = 88°, ro = 110 At/, aout ~ -5, a/„ = 2, g = 
0.50, corresponding to the Ks band. Red lines respectively blue 
lines stand for the NE side respectively the SW side. Dotted lines 
show the disk intensity dispersion at 1 cr within the slit window 
of 0.25". Dashed lines correspond to the fitted power law profile 
next to which the slopes are indicated. 

7. Limit of detection to point sources 

Our data feature the highest contrast ever reported for HD 32297 
in the Ks band near ~1" hence suitable to derive limit of detec- 
tion of point sources in the planet/brown dwarf mass regime. In 
this particular case, it is relevant to measure the contrast along 
the midplane as the orbits of hypothetical planets will be likely 
aligned with the disk plane. For each image in Fig. 1 the disk 
is isolated in a 0.3" slit window centered on the midplane. The 
standard deviation is calculated in the pixels that are vertical 
to the midplane at each radii (we check that the result do not 
strongly depends on that slit width). The photometric biases in- 
troduced by the self- subtraction of the three ADI algorithms 
were calibrated with nine fake planets (copy of the PSF) injected 
into the data cube and located at specific locations (from 0.25" 
to 5" and 90° from the midplane to avoid confusion with the 
disk). An attenuation profile is interpolated on these flux losses 
assuming that it varies monotonically with radius (not necessary 
true for LOCI) and used to correct for photometric biases and 
to provide calibrated contrast curves. The 5cr contrast curves are 
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Fig. 10. Detection limit to point sources at 5cr for cADI ( ma- 
genta), rADI (red) and LOCI (blue). The expected intensity of 
some planetary objects from the BT-SETTL model are shown 
for an age of 30 Myr (dashed lines). The dotted lines corre- 
spond to the mass/separation range of planets able to sculpt the 
disk and assuming that the inner edge of the ring is located at 90, 
100 and llOAU. 

presented in Fig. 10. As it is measured in the disk midplane in- 
stead of the background noise, the limit of detection is similar 
for the three algorithms. We plot also the level of contrast for 
three planet masses: 6, 10 and 20 Mjupiter assuming an age of 30 
Myr (Kalas 2005) and considering the BT-SETTL models from 
Allard et al. (201 1). However, we note that the age of HD 32297 
is not very well determined. We can definitely reject the pres- 
ence of brown dwarfs at separations larger ~50 AU and closer 
than -600 AU (FOV limitation). Due to the distance of the star 
(1 13 pc) it is impossible to put any constraints on planetary mass 
objects at closer separations. At 1 10 AU, the distance of the inner 
ring, the detection limit reaches -lOMjupiter- Assuming the inner 
cavity is sculpted by a planet and following the approach pre- 
sented in Lagrange et al. (2012b) we can estimate the position 
of that hypothetical planet. According to Wisdom (1980) and 
Duncan et al. (1989, Eq. 23) the distance of the planet to the rim 
of the ring (Sa) is related to the planet/star mass ratio (Mp/Ms) 
with the relation: 

Sa/a= \3x{MplMsf'^ (3) 

where a is the semi-major axis. Soa+6a is the distance of the in- 
ner rim which is lower than tq. The same formulation is used by 
Chiang et al. (2009) in the case of Fomalhaut b. We calculated 
Mp for several assumptions on a + 6a (90, 100, llOAU) and 
taking Ms = 2.09Mq. We interpolated these masses on the BT- 
SETTL model, hence converted into absolute magnitudes and 
into contrasts. Smaller values of a -\- 6 a would lead to unrealistic 
masses in the stellar domain. The dotted lines in Fig. 10 show 
the combination of planet masses and separations that are able 
to produce an inner ring at llOAU. The comparison between 
dynamical constraints (dotted lines) and detection limits (solid 
lines) indicates a possible range of separations and masses for 
the planet: a > 65 ~ 85 AU and Mp <12 Mjupiter- We stress 
that this estimation is only valid for a planet at the quadrature 
(projected separation = physical separation). However, given the 
steepness of the dynamical constraints, if a planet produces the 
ring it is necessarily close to the rim or it would have been mas- 
sive and hence detected. 



8. Summary & conclusions 

We present high contrast images of the debris disk around 
HD 32297 obtained with NACO at the VLT in the near IR H 
and Ks bands. To improve the attenuation of the starlight and the 
speckle rejection we combined phase mask coronagraphy and 
ADI algorithms. The disk nearly edge-on is detected in the two 
filters at radial distances from 0.65" to 2" (74 to 226 AU at 
the distance of 113pc) owing to various ADI algorithms. ADI 
is a very efficient method to achieve detection but induces some 
photometric biases that have to be corrected for in the data anal- 
ysis for a proper interpretation. We observed a deviation from 
the midplane near 0.8" in the NE side that is clearly not an ADI 
bias. We focused our eff'ort on the morphological modeling of the 
disk using a numerical code GRaTer, which produces synthetic 
scattered light images of debris disks. Following previous ob- 
servations, we considered a disk containing a planetesimal belt 
with dust distribution inward and outward. We generated a 5- 
parameter grid of models and applied several sorts of minimiza- 
tions to constrain the disk morphology. 

As a first result, we found in Ks and H bands a higher in- 
clination than previously published for l.l/im (88° instead of 
80°). In our images, the vertical extent of the disk together with 
the disk parameters found from modeling are definitely not con- 
sistent with such a disk inclination. However, we suspect that 
the isophotal ellipse fitting used in former data (Schneider et al. 
2005) can be biased by the strong interaction with the inter- 
stellar medium as revealed by Debes et al. (2009), especially 
when considering the lower angular resolution and conversely 
the higher sensitivity to faint structures of NICMOS/HST as op- 
posed to NACO. In addition, we did not found any significant 
NE/SW brightness asymmetry like in Schneider et al. (2005) and 
Mawet et al. (2009) although we are also observing at near IR 
wavelengths too. 

Furthermore, the minimization of the disk parameters al- 
lowed to constrain the radius of the planetesimal belt at 1 10 AU. 
A NE/SW asymmetry is suspected, possibly due to an off'set of 
the planetesimal belt with respect to the star. We showed that the 
deviation from the midplane is likely the result of anisotropic 
scattering. One edge of the disk is therefore brighter. Whether 
it is linked to the properties of the dust or to the geometry 
(the hypothetical off'set of the belt) and hence to a gravitational 
perturber is yet undetermined. It is remarkable that among the 
known debris disks, three of them surrounding A-type stars, 
namely, J3 Pic, Fomalhaut, and HD 32297 share the same posi- 
tion of the planetesimal belt near 100 AU as a possible indication 
of a common architecture. 

Intensity related parameters like the brightness slopes of the 
inner and outer parts are more difficult to constrain and the val- 
ues we determined should be taken with caution as the uncer- 
tainty can be large. However, it gives an order of magnitude 
indicating that the outer parts have a power law index close to 
that expected for dust grains expelled from the system by radia- 
tion pressure. As for the inner parts, there is a significant diff'er- 
ence between Ks and H band filters, but overall, the steepness is 
smaller than some other ring-like debris disks as already noted 
by Moerchen et al. (2007). 

Armed with a detailed morphological characterization it is 
now possible to start a dynamical study to understand the prop- 
erties of the debris disk around HD 32297 in regards of planetary 
formation or gravitational perturbation as it was performed in 
other cases (Wyatt 2005; Chiang et al. 2009). Finally, the obser- 
vations of debris disks from ground-based telescopes in the near 
IR are achieving comparable or even better angular resolution 
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than from HST, owing to adaptive optics. A similar statement 
can be made about contrast owing to a combination of high con- 
trast imaging techniques (coronagraphy, ADI). The next genera- 
tion of extreme adaptive optics instruments, the so-called planet 
finders (Beuzit et al. 2008; Macintosh et al. 2008; Hodapp et al. 
2008), will not only revolutionize the search and characteriza- 
tion of young long-period planets but also the field of debris disk 
science pending the arrival of the James Webb Space Telescope 
(Gardner et al. 2006). 
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